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We analyze the nonequilibrium Kondo model at finite voltage and temperature by using a new 
formulation of the real-time renormalization group method with the Laplace variable as the flow 
parameter. We evaluate the energy-dependent spin relaxation rate and nonlinear conductance, 
and derive an approximate form for the universal line shape for the latter in the whole crossover 
regime from weak to strong coupling. The results are shown to agree well with exact methods 
in equilibrium, Fermi-liquid theory, weak-coupling expansions, and recent experiments. For the 
transient spin dynamics we find a universal exponential decay in the long-time limit along with a 
truncation-dependent pre-exponential power law. For multichannel models a pure power-law decay 
typical for non-Fermi-liquid behaviour is predicted. 
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The solution of the nonequilibrium Kondo model is 
a basic unsolved problem in condensed matter physics. 
In its elementary version, the model consists of a spin 
^ coupled antiferromagnetically to two fermionic metal- 
lic reservoirs. It reveals the Kondo effect of complete 
spin screening at low energies and plays an important 
role in bulk systems with magnetic impurities [l[ . In the 
Coulomb blockade regime of single-level quantum dots, 
the Kondo model can be realized via its relation to the 
symmetric Anderson impurity model |3- In this case the 



effective exchange coupling is given by J 



2A 

TlU- 



where 



A is the level broadening and U ^ D is the charg- 
ing energy, defining the band width 2D. This allows 
for a controlled experimental realization of the Kondo 
model in nonequilibrium, where the Kondo effect has 
been observed via a zero-bias anomaly Q. In equi- 
librium, the Kondo model has been solved by a Bethe 
ansatz y and conformal field theory (CFT) [5|- Using 
the numerical renormalization group (NRG) method [gi- 
ll] and a Bethe ansatz [§], the temperature dependent 
linear conductance G{T) has been shown to be a uni- 
versal function of T/Tk, where Tk De-^/^^-^) is the 
Kondo temperature. Analytical results are available for 
low (T <C Tk) and high (T Tk) temperatures from 
Fermi-liquid theory and poor man scaling methods [lol - 
[12I 1 . A challenging issue is the nonequilibrium properties 
at T = 0, such as the determination of the nonlinear 
conductance G{V) as a function of bias voltage V and 
the transient spin dynamics {S}{t). The weak-coupling 
regime max{F, 1/t} ^ Tk has been solved by improved 
poor man scaling methods using phenomenological decay 
13 flow equations O, [i3|j ^'^id the real-time 
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renormalization group (RTRG) method 
ever, the analysis of universal properties in the strong 
coupling regime V, 1/t < Tk remains an open problem. 

In this Letter we propose an approximate solution of 
this problem by using a new formulation of the RTRG 
method where the Laplace variable E is used as the flow 



parameter (called "i?-flow" in the following) . Within this 
scheme the relaxation rates occurring in each renormal- 
ization group (RG) step are the full physical ones, and 
their energy dependence appears to be crucial for the de- 
scription of the crossover from weak to strong coupling. 
Applying the £'-flow to the nonequilibrium Kondo model, 
we calculate the nonlinear conductance and observe that 
both the T and V dependence agree well with recent ex- 
periments [2l[. We show analytically that in the scaling 
limit {D — >■ 00, J — >■ 0, Tk = const) both the low and 
high energy properties are consistent with Fermi-liquid 
theory [ll| and weak-coupling expansions, respectively. 
This distinguishes our approach from other studies of 
the nonequilibrium Anderson impurity model (see, e.g., 
Refs. and [2^ for reviews), where the effective spin ex- 
change coupling J cannot be chosen arbitrarily small to 
achieve universality in the scaling limit. Also approx- 
imate calculations of G{V) via the Bethe ansatz eigen- 
states [§| cannot cover the crossover to the weak coupling 
regime. We find that our G{T) curve agrees within a few 
percent with the NRG one and that the nonlinear con- 
ductance G(y) is almost independent of the truncation 
order. 

Finally we address the transient spin dynamics after 
an initial coupling of the spin to the bath. In Ref. [20| it 
has been predicted in the weak coupling regime that the 
long-time dynamics is always exponential, accompanied 
by pre-exponential power laws in higher orders. Scal- 
ing behavior has been observed within time-dependent 
NRG calculations and exponential behavior was pre- 
dicted at the Thoulouse point [2^ and in the frame- 
work of integrable quantum field theories (26j . How- 
ever, due to the presence of branch-cut contributions 
the transient spin dynamics is still under debate (see 
Ref. 27| for a discussion within the spin boson model). 
At T = y = we find a universal exponential decay in 
the long-time limit t 3> 1/Tk with a rate F* ~ Tk, and, 
in addition, a pre-exponential power-law behavior with 



2 



a truncation-dependent exponent. In contrast, for multi- 
channel Kondo models, we predict a pure power-law de- 
cay in the long time limit as a result of non-Fermi-liquid 
behaviour at low energies d, [2^ . 

Model. — We consider a single spin ^ with spin S_, which 
is coupled by an isotropic exchange H^x — JoS_ ■ s to the 
spins s = i Eaa'^^'fefc' aLfe ^<7<7' ««'a'fc' of two reservoirs 
labeled hy a = L,R. a =i,t denotes the spin, k is the 
state index and a are the Pauli matrices. The reservoirs 



are described by Hr, 



with finite 



band width and are kept on chemical potentials Hl/ji — 
±V/2 (we use units e = ?i = fcs = 1). _ 

Liouvillian approach. — Following Ref. lisl . we describe 
the dynamics of the reduced density matrix p{t) of 
the local spin by the effective von Neumann equation 
p{t) = —iJ^dt'L{t — t')p{t'), which reads in Laplace 
space p{E) = i/{E - L{E))pQ At the initial time 

t — the system is assumed to factorize into an arbi- 
trary local part and a grand canonical equilibrium 
part p°^p°^ for the reservoirs. L{E) is an effective Liou- 
villian, which results from integrating out the reservoir 
degrees of freedom. Using spin conservation and rota- 
tional invariance for the isotropic Kondo model we find 
that the Liouvillian can be parameterized by the spin re- 
laxation rate r(i?) via the nonvanishing matrix elements 
= -^ss'T{E) and L{E)s,-s,s,-s = ~ir{E), 
where s =tii= i- The Liouvillian has a zero eigenvalue 
corresponding to the stationary state and three degen- 
erate eigenvalues —iT{E) describing the spin dynamics 
via 



(S){E) 



E + iT{E) 



(5>(t = 0) 



(1) 



from which the real-time dynamics can be obtained 
via the inverse Laplace transform. Analogously, one 
can also derive an equation for the current Ia{E) = 
—iTvYia{E)p{E) fiowing into the reservoir a, where 
Yia{E)ss,s' s' = i^aiE) are the nonvanishing matrix 
elements of the current kernel. Defining G(E) = 
-K-^V i^(Ei)^ the stationary conductance in units of Go = 
2eV/i is given by G = G(0). 

RG formalism. — To calculate T{E) and G{E) we use 
a different formulation of the RTRG method than the 



one proposed in Ref. Il8l . Instead of introducing an arti- 
ficial cutoff A on the imaginary frequency axis, we use 
the physical Laplace variable E itself as the flow pa- 
rameter. This means that we differentiate the full di- 
agrammatic series for the effective kernels and vertices 
w.r.t. E and obtain a self-consistent set of equations in 
terms of effective quantities by an appropriate resumma- 
tion of all diagrams (see details in the Supplementary 
Material fs^). It turns out that the scaling limit is well- 
defined if two iJ-derivatives are taken for the kernels, 
whereas a single derivative is sufficient for the vertices. 
The exact RG equations are finally truncated systemati- 
cally by considering all diagrams up to a certain order in 



the effective vertices but keeping the full propagator be- 
tween the vertices. The reliability of our approach in the 
strong coupling regime is tested by comparing the results 
in second and third order truncation. These truncation 
schemes are essentially different from conventional 1- and 
2-loop treatments, since we consider the full propagator 
between the vertices and since we use the same trunca- 
tion scheme for the rate T{E) including its full energy 
dependence. We illustrate this by showing the RG equa- 
tions for T = F = (the full equations at finite T, V can 
be found in the Supplementary Material |30(]). Neglecting 
all terms of 0{J'^), we obtain with E = iA 



dlr 



where Z 



4.P 

' A + T 



dKJ = - 



2J2(1 + ZJ) 

A + r 



(2) 



denotes the Z-factor. With A = A 



the RG equation for the conductance reads c^aG = 
— ^^JjK, which is coupled to two other vertices with the 
RG equations d\Ji — — 2J/J/A and d^K = — 47^J(1 — 
ZJ)/X. With J = ZJ wc can write Eq. dU in the 
form dxJ = — 2J(1 — J)/A which gives the invariant 

Tk = (A -I- r)^ YZj'e^^j defining the Kondo temper- 
ature. In contrast, the standard 2-loop Kondo tempera- 
ture is obtained by neglecting the relaxation rate F lead- 
ing to the poor man scaling solution Jp(A), defined by 

Tji = Ay^Y^-e '^''p . Thus our approach provides a mi- 
croscopic foundation for the physically intuitive picture 
of the relaxation rate cutting off the RG fiow, including 
its full energy dependence. As an important consequence 
wc find in our approach that the 2-loop fixed point at 
A = is shifted to the negative value A = — F* < 0. 
As shown below this is fundamentally related to the ex- 
ponential relaxation of the spin dynamics. Moreover, a 
technical advantage is that the point E = iA = 0, defin- 
ing the stationary value of all physical quantities, is well 
separated from the fixed point. This allows for a well- 
defined expansion around this point leading to the cor- 
rect Fermi- liquid relations (see below). 

The initial conditions at high energies Ao ^ Tk are 
given by J = J07 ^ = 1, ^ = 2Jq, J/ = Jq and 
G = ^^Jq ^M- the scaling limit, this leads to the 
solution Z = {1- Jf and G = ^ For F a special 
procedure is needed since it is proportional to A in the 
weak coupling regime F = 2A( Jp -I- 0( J^)). This leads to 
a truncation order independent instability for its value 
in the regime A ^ Tk. We solve this problem by com- 
plementing our approach with the exactly known result 
of unitary conductance G{E = 0) = ^J'^{E = 0) = 1, 
which is a consequence of the Friedel sum rule [l|. This 
fixes the values J = J(0) 



^ - 0.37, f ^ F(0) 



i-j^e27 Tk « 5.11 Tk, and Z = [l - Jf k 0.40, and 
the RG equations © can be solved numerically starting 
from A = 0. The T = V = Q results for all quanti- 
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FIG. 1. (Color online) Main panel: G{T) in second (red 
dashed line) and third (blue full line) order truncation as 
function of T/T^, defined by G{T = T|-) = 1/2, compared 
to the NRG results (green diamonds). Inset: (Grtrg — 
Gnrg)/Gnrg as function of T/Tj^ in second (lower line) and 
third (upper line) order truncation. 



ties at large Aq 3> Tk are then used as initial condition 
to solve the RG equations at finite T, V using an initial 
value Ao>T,l/,Tif. 

Finite T. — As a benchmark for the reliability of our 
approach we first compare the linear conductance G{T) 
with the NRG result. When plotted in units of T^, de- 
fined by G(T = T^) = 1/2, we see in Fig.[T]that our re- 
sults in second and third order truncation fall almost on 
top of each other and agree very well with the NRG result 
. The low- and high-temperature regimes can be stud- 
ied analytically. For T > TJ^, we expand in Jp(T) ^ 1 
and get 



FIG. 2. (Color online) G{V) in second (red dashed line) and 
third (blue full line) order truncation as function of T /T^ , 
defined by GiV = Tj") = 1/2, compared to Lorentzian fit 
(green squares) and high-!/ result (Q (dotted-dashed line). 



of a term in 0( J^). For small V < T^*, we get G{V) « 
1 — cy(^)^ with cv/ct = 3/(27r^). This coincides with 
the exact ratio known from Fermi-liquid theory [ill,!!!. 
The result is fulfilled order by order by expanding in J. 
The numerical solution does not show the exact ratio 
since inconsistent terms of 0{J^) are involved there, but 
in third order truncation the error is already below 10%. 
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In Fig. [2] we show a comparison with a Lorentzian 
fit function G{V) = 1/[1 + cv{V/f)% as proposed in 
Ref. 9] from Bethe ansatz calculations. It turns out to 
be a good fit only for voltages significantly below T^*, 
but becomes inconsistent with the weak-coupling solu- 
tion when applied in the whole crossover regime. The al- 



G(T) « —Jp{T) {^1 - Jp{T)A\n--^ + 0{J^) j , (3) ternativc trial function G{V) = [1 + {V/T'j^{x)y 

where 7 = 0.577 ... is Euler's constant. Besides the stan- 
dard weak-coupling result, the second term is a third or- 
der correction which suppresses the conductance. For 
T <C Tk, we expand in J and obtain in third order trun- 
cation the Fermi-liquid result G{T) w 1 — ct(^)^ with 

CT = TT^ '^l[t^f^3^ + 0{P). For the absolute value of the 
curvature w.r.t. Tj^ we obtain in third order truncation 
CT(Tff/r)^ « 4.8, which is much better than the second 
order result but still not fully reliable due to the inac- 
curacy in the ratio T^/f (the NRG result is about 10% 
larger 31]). 

Finite V. — Next we consider T — but finite V. Fig- 
ure [2] shows again that our second and third order trun- 
cation results for G{V) agree very well with each other 
when plotted in units of T^* defined by G{V = T^) = 
1/2. For large V > T^* we expand in Jp{V) < 1 and 
obtain 



3^' 



G{V)^—Jl{V) + 0{jf,) 



(4) 



The first term is the standard weak-coupling result, but 
in addition our third order analysis proves the absence 



with 

T'j^{x) = T^*(if^^rLL)i/2 and x = V/T*^*, which is a 
good fit for G(T) at s = 0.22, & = and T^* ^ 
6], turns out to describe the crossover regime of G{V) 
quite well at the significantly different values s — 0.32, 
h = 0.05 and s' = 1.26. This shows that the univer- 
sal line shapes of G(y) and G(T) are quite different and 
cannot be rescaled on top of each other, as shown in 
Fig. [3] when plotted in terms of the same unit T^. A 
numerical analysis gives T^*/T^ ~ 1.8 and we find that, 
in third order truncation, G{V) lies always above G(T) 
for V = T , which can already be anticipated from our 
analytical results at low and high energies. In Fig. [3] 
we also show that recent experiments appear to be well 
captured by the present theory. The experiments are 
performed on an InAs nanowire quantum dot [2l| , which 
can be represented as a single-level dot with the effec- 
tive spin exchange J ^ ^ ^ 0.28. The experimen- 
tal data for G{T) and G{V) fit our results quite well 
in the regime T,V using the same rescaling fac- 

tor for the V and T dependence. Since ~ 0.1 J7 
in this experiment, universal features from spin fluctu- 
ations can only be expected below T^, while above Tj- 
charge fluctuations will set in. As proposed in Ref. [33|, 
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FIG. 3. (Color online) Comparison between the RTRG results 
and experimental data [IJ for G{T) and G{V). Both T and 
V are expressed in units defined by G{T = T^) = 1/2. 
An asymmetry factor AxrXl ~ 0.8 is introduced to match the 
experimental G(0) with AxrXl- 

the universal relation G{V = T^) w | allows for an el- 
egant experimental determination of the important en- 
ergy scale T^. The experimental conductance is not 
unitary due to uncontrolled asymmetric couplings to the 
two leads. We have checked that asymmetric couplings 

Jaa> = I^XaXa'J with < Xq, < 1 and Xh + = 1 

only give rise to the well-known asymmetry factor Axlxr 
for the conductance in the weak and strong coupling 
regime JJ), 32] and change the normalized conductance 
G{V)/{'ixLXR) in the intermediate voltage regime by less 
than 1% for the experimental value AxlXr ^ 0.8. Fur- 
thermore, we checked that the finite experimental tem- 
perature T - lOmK - 0.025 does not change G{V) 
significantly. 

Time evolution. — At T = V" = and for t ^ 1/Tk 
we find from Eqs. ([T]) and ^ that the spin dynamics is 
given by {S){t) - j^e"^'*. The resolvent l/[E + iT{E)] 
has a branch cut on the negative imaginary axis with 
a branching pole at _E = — iF*, which is identical to the 
fixed point J = J* of the RG flow and determines the rate 
of the exponential decay. We obtain the nearly constant 
ratio T* /Tk ~ 1.1 independent of the truncation order. 
To determine the exponent g we expand around the fixed 
point and set E = iA with |A -h r*| < F*. We obtain 
^ = 4Z JVA w aZ/X with a = 4J*^. This gives Z - A" 
and, by using d^X = Z^^ ^ A^", we get A ^ (A -|- 
p*-)i/(i+a)^ leading to g = = after the inverse 

Laplace transform. However, the fixed point value J* 
depends on the truncation order. In second (third) order 
truncation we get J* = 1 and g = | {J* = 5 and g — 
i). As a result we obtain a universal exponential decay 
but the prc-exponential power law cannot be determined 
unanibigiously. 

In contrast, for iV-channel Kondo models with iV > 1, 
where ^ = ANZJ^/X and ^ = -2J2(1 - 7VJ)/A, the 
non- Fermi- liquid behavior at low energies [^, [2^ requires 



the fixed point to be at = 0. This gives rise to a 
pure power-law decay for the spin dynamics in the long- 
time limit. The exponent follows from a — ANJ*'^ and 
J* = to be g = ^^^^'j,'2 = However, the precise 

value of the exponent is only reliable for N ^ 1, where 
the /3-function can be systematically truncated [34|. 

Conclusions. — We have derived an approximate form 
for the universal line shape of the zero-bias anomaly for 
the isotropic Kondo model, which successfully describes 
the recent experiments [U [s^. The reliability of our 
approach has been checked by comparing different trun- 
cation orders and by reproducing the correct asymptotic 
values at large and small voltages. In the long-time 
limit we found a universal exponential decay for the spin 
dynamics in the 1-channel case and a power-law decay 
for multichannel models. Our approximate results could 
serve as a reference for future more refined numerical and 
analytical approaches. In particular, the precise value 
of the exponential decay rate and the form of the pre- 
cxponential modulation is an open issue. The new for- 
mulation of the RTRG method using the Laplace variable 
as the fiow parameter is a promising analytical tool for 
solving other nonequilibrium problems of quantum im- 
purity physics and dissipative quantum mechanics in the 
strong coupling regime. 
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The following appendices comprise the supplementary material to Ref. [l|. 



I. THE S-FLOW SCHEME OF RTRG 



In Ref. 12 a real-time renormalization group method (RTRG) has been developed in Liouville space, which can be 
used to evaluate the effective Liouvillian L{E) of a local quantum system coupled to several fermionic reservoirs. 
Within this method, a Matsubara cutoff was defined by cutting off the Matsubara poles of the Fermi functions. A 
systematic weak-coupling expansion including cutoff scales from relaxation and decoherence rates was developed, 
which was used to calculate the nonequilibrium properties of the Kondo model^ and the interacting resonant level 
model (IRLM)* in the weak coupling regime. In the strong coupling regime and in the presence of spin, orbital or 
potential fiuctuations, a problem within the RTRG method is the generation of spurious terms linear in the physical 
band width 2D of the reservoirs, originating from the fact that the cutoff used in the RG is not identical to the physical 
band width. In particular, it turns out that the RG-flow is unstable against a small variation of these terms. In the 
weak-coupling regime it has been showi*2ii^ that the problem can be resolved by means of an appropriate subtraction 
scheme. The same was demonstrated for the strong coupling regime of the IRLM**. However, in the generic strong 
coupling case it is not clear whether the subtraction scheme can be uniquely defined. 

Here we propose an alternative fiow scheme where this problem is avoided by solving the RG-equations starting 
at low energies. It uses the Laplace variable E as fiow parameter and for this reason is called i?-fiow scheme in the 
following. Another advantage of the £'-flow scheme consists in the fact that the Laplace variable is a physical energy 
scale from which the time evolution can be directly evaluated. The effective Liouvillian and the vertices occurring in 
the RG equations are the full physical ones. Furthermore, the method allows the consideration of all diagrams up to 
a certain order without adhering to approximations typically used in the weak-coupling regime. 

We show here the derivation of the £'-flow RG equations in third order truncation for the case where only fermionic 
2-point vertices Gi2(£', wi, W2) are present and describe either spin, orbital or potential fiuctuations (higher-order 
vertices are only generated in fourth order). Moreover, we assume that the bare vertices have no explicit frequency- 
dependence and we take a flat d.o.s. in the reservoirs with band width 2D. The index 1 = rjaa characterizes 
the reservoir field operator: rj — ± distinguishes creation/annihilation operators, a is the reservoir index, and a is 
the channel index (i.e. the spin for the Kondo model). The frequency arguments uii = rjiUJi contain the reservoir 
energy measured relative to the chemical potential fj,a of reservoir a. The Keldysh index p can also be included 
in the index 1 but is disregarded here since, as will be shown below, only the average over the Keldysh indices 
remains in the final RG equations. Including the Keldysh indices the contraction between two vertices is given by 
7^^, (cj, oj') = Sii>5{ui + Cd')^P (w), where 7^ (cD) = p'f{p'id) j^P^^i , 1 = "qaa, 1' = —ij'a'a', and /(w) = e--^/^+i denotes 
the Fermi function at temperature T. 

To achieve convergence of the frequency integrals in the scaling limit D ^ 00 one needs equations for the second 
derivative of L(E) and the first derivative of Gi2{E, uji , 0J2) w.r.t. E. Using the diagrammatic representation described 
in Ref. i2i, the i?-derivative can only act on the propagators H(i?) = occurring between the bare vertices (we 

always resum all self-energy insertions such that the full effective Liouvillian appears in the propagator). Resumming 
all diagrams right and left to the propagator in which the £'-derivative is taken, we obtain the following two equations 
up to 0(G3) 



— 612(^^,^1,^2) 



5)6 



(102; 



0{G^) 



(2) 



Here the slash indicates the i?-derivative and a double-circle represents the full effective 2-point vertex (containing 
all connected diagrams with two free reservoir lines) . Symmetry factors ^ arising either from the diagrammatic rules 
(when two vertices are connected by n equivalent lines) or from the ^^-derivatives -^d^ are explicitly quoted. Since the 
propagators Hi...„ — H(£'i...„ + uJi...n) between the vertices contain the variables i?i...n — E p,i...n = E + p,i + . . . p,n 
and wi...n ^ LJi + . . . LJn with [li — r]i^a. and uii = r]iUJi, the i?-derivative can be written as frequency-derivative 
and we can apply partial integration to calculate the frequency integrals. E.g., for the first term on the r.h.s. of ^ we 
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obtain (we have permuted the two indices of the vertices by using antisymmetry Gi2{E,uji,uj2) = ^G2i(£',a;2, wi)) 



(1)6 ^ 6o i6-^^-4(>^ — (5 (5<) 



(3) 



12 12 12 

Here the cross indicates the frequency derivative w.r.t. uj of either the corresponding reservoir contraction or the 
frequency argument of the vertices. The first hue is exact and follows from partial integration. For the derivation of 
the second line we have used 



^ = ic^^Z^ + O(G^) , ^ = =J^^ + 0(G3) , (4) 

which Vollows analogously to y using the fact that, in the original diagrammatic series, a frequency associated 
with an external line occurs only in those propagators which lie below the external line but not in other propagators 
and not in the bare vertices. Note that these relations can only be applied if it is specified whether the external line 
involving the frequency derivative is directed either towards the left or towards the right. 

Analogously, we can treat the first term on the r.h.s. of ([T]) by two partial integrations and using the fact that 

-Gi2{E, Qi, Q2) - 0{G^). We ohtam 



(5) 



where, in the second step, we have again used (|3]). 

Inserting ([5]) and Q in ([Ij and ([2|), respectively, we see that many diagrams in 0{G^) cancel each other and we 
cast the final RG-equations to a very compact and generic form 



L{E) 



0{G^), 



(6) 



d_ 

dE 



Gi2iE,UJi,UJ2) 



-(lo2) 



00 

1 2 



3) + OiG% (7) 



1 2 

The RG equations are constructed in such a way that all frequency integrals on the r.h.s. of the RG equations 
are well-defined in the limit Z? — > 00 even when the frequency dependence of the vertices is neglected. Performing 
this limit, each reservoir contraction contains (ui) — p'f{p'ui). As a consequence, the frequency integral for a 
contraction where the derivative is taken (diagrammatically indicated by the cross) involves (w) = /'(w), which 

is independent of the Keldysh indices. The frequency integral of the contraction without a derivative in the second 



term on the r.h.s. of (O involves 7^ (w) — + y, where = ,f{u 



is the antisymmetric part of the Fermi 



function. Since two propagators are involved in this frequency integral, the integration contour can be closed in the 
upper half of the complex plane and the symmetric part ^ vanishes identically since the resolvents and the vertices 
are analytic functions in the upper half planed. Therefore, an explicit dependency on the Keldysh indices disappears, 
and it is sufficient to consider the vertices Gi2{E, 0)1,012) — X^pipa ^12^^ '^i' ^^2) averaged over the Keldysh indices. 
We note that the symmetric part of the Fermi function does not influence the RG equations but enters into the initial 
conditions at high energies^. 

To evaluate the frequency integrals on the r.h.s. of the RG equations (|6]) and (O explicitly, we need a consistent 
approximation for the frequency dependence of the vertices and the Liouvillian. For the vertices, this can be achieved 
by formally integrating the first two terms on the r.h.s. of ([7]) yielding 



I — o 



-o- 



i — t: 



0{G^) 



(8) 



1 2 2 1 

where Gi2{E) = Gi2{E,0,0) is indicated by filled double dots. The open circle at a contraction indicates that 
the propagator corresponding to the vertical cut at the position of that circle has to be replaced by the difference 
n(i?i...n + + w) — n(i?i...„ + 0)1. ..n) of the two propagators with and without the frequency ui corresponding to 
the contraction with the circle. The identity refers to the case where the two external lines are directed to the 
right, similar equations can be written down for the other cases. We note that the sign factor in the second term on 
the r.h.s. accounts explicitly for the crossing of the two external lines relative to the first term. Therefore, when using 
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(IHl) in a certain diagram the sign factor must not be written explicitly since it is automatically accounted for in the 
diagrammatic rules. Inserting ©in ((5]) and (0 we obtain 



L{E) = 



(9) 



dE 



(i«2; 









1 — X — 






(>i — . 


H 


. i<. 



1 2 



(10) 



Finally, to calculate the frequency integrals in these equations, one needs a consistent approximation for the propagator 

n(i? + (D), where u) is an integration variable. Using ■^^L{E) ~ 0{^G^) cf. (|9]), we approximate L{E + w) « 

a 



L{E) 



dE 



L{E)u. This leads to the following approximation for the propagator 



n{E + Ld) 



X(E) 



Z{E) 



(11) 



where Z{E) 



is the Z-factor operator and x(^) = ^{E){E — L{E)) is an operator defining the distance 



to resonance positions. We note that the eigenvalue zero of the Liouvillian L{E) (or of Z{E)L(E)) can be omitted 
in pT|) since the projector Pq on the corresponding eigenvector fulfils PqL{E) = -P0G12 = and PqZ{E) = 1—. Using 
PT|) in ^ and pUl) , all frequency integrals can be calculated analytically, which will be done explicitly for the Kondo 
model in the next section. 

Equations ([9]) and pO)) along with ([TT|) are the final RG equations for a determination of the effective Liouvillian 
L{E) from which the reduced density matrix of the local system in Laplace space follows via p{E) = Po, where 

po — p(t — 0) is the initial condition. To find the average of the current la — —iTTT,a{E)p(E) flowing into reservoir 
a, one needs the current kernel 'Sa{E) in Laplace space. As shown in Ref. 2 it can be determined analogously to L{E) 
by replacing the first vertex from the left in all diagrams of ^ and (|10l) by the current vertex Ii2{E). To calculate 
the differential conductance one needs the variation SIq, for an infinitesimal variation 6pa of the chemical potentials of 
the reservoirs. Within the _B-flow scheme the RG equation for ■^SL{E) (or, equivalently, for -^ST^aiE) by replacing 
the first vertex by the current vertex) can be straightforwardly established by applying the variation to the original 
diagrammatic series and using SIIi,,,n = 6fLi,,,n-^^i...n + ni...„(5Li...„ni...„, where = (5L)(i?i...„ + a'i...„). 

After resummation we obtain by analogy with ([T]) up to 0{G'^) 



1 2 



(12) 
12 3 

where SL ^ 0{SpG) is represented by ©— ^^'l '^^'^ ^'^ calculated from the first (lowest order) term in the r.h.s. 

of ()12p and subsequently inserted in the second and third term in the r.h.s. of ()12p . Applying two partial integrations 
to the first term on the r.h.s. of (fT2|) by analogy with ([5]), we obtain 



_d_ 

dE 



SHE) = - l^^l 



0{G^) 



(13) 



Using ([8|) we get the final RG equation for the variation of the kernel 




5pi2 ii^^ir^"^ 



O(G^) 



(14) 



13 1 i 

We note that at zero temperature all diagrams in ^ , pU)) and containing a contraction with a circle and a cross 
vanish, since the cross restricts the frequency to zero value. 
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II. RG EQUATIONS FOR THE KONDO MODEL 

We now apply the RG-equations (|9]), (fTO)) and (|T4l) to the isotropic spin-i Kondo model at zero magnetic field. 
Using spin conservation and rotational invariance it has been shown in Ref. [H how the kernels and the vertices 
can be decomposed into convenient basis superoperators. One obtains L{E) = —iT{E)L°' , 'Ej{E) = iT^{E)L^, 

G+a.,-aV'(£;) = E^=a,2,3GL'(^^)^L', ^nd I^^.^^,,, (£;) - E;,=M ^-'^ (^)^-" ^here L^;,^ = L^H„„, and 
L^g., — V^crcr' fo'" * = Ij 2, 3 (ct are the Pauli matrices). The superoperators L°'^ and L^'^''^ can be expressed in terms of 
the spin superoperators L+ = S- and = -{-S) via = \+L^L~, = j-L^L^^L}-^ = ^{L+ -L~ T'^iL'^ AL~), 
and = —^(L^ + LT)- The closed algebra of these superoperators needed to evaluate the RG-equations can be 
found in Ref. The algebra is applied after the summations over the creation/annihilation and spin indices rj and 
a have been performed. Furthermore, by using dill) , the intermediate propagators are replaced by the c- numbers 
I1{E + w) w u>lf(E) , where Z{E) = ^^^^^^^ and x{E) = Z{E){E + iT{E)). This follows from the fact that the zero 

eigenvalue of the Liouvillian can be omitted, and the three nonzero eigenvalues are degenerate and given by —iT[E). 

A detailed analysis shows that the vertex component /'''''' is generated but does not influence the current kernel. 
Furthermore, with ^ J, it turns out that the vertex component G" is generated in O(J^) and influences the 
Liouvillian only in O(J^). Therefore it can be omitted as well and we are left with a set of closed equations for the 
vertex components G^, and 7'''^^, together with the rates F and F^. For applied voltages the energy arguments 
are shifted by multiples of the chemical potentials of the reservoirs and all components couple to each other. After 
the summations over the indices rj and a have been performed, we use the shorthand notation 1 = a for the reservoir 
index. Furthermore, we define /ii2 = jiai — fJ-a2 1 ^12 — E + /ii2 and -B1234 ^ E + /ii2 + /i34. For convenience we define 
Ji2{E) = -Gi2{E), Ki2{E) = -i|G?2(^^2i) and I^2iE) = -41^2^^). From hermiticity of the original Hamiltonian, 
we obtain the symmetry properties Ji2{E)* = J2i(--E*), Ki2{E)* = K2i{-E*) and /72(^)* = --^2i(--^*) for t^e 
vertex components, and T{E)* = T{—E*) and T^(E)* — T^(~E*) for the rates^. The initial conditions at high 
energies arising from the symmetric part of the Fermi function are calculated in Ref. !2 as 

a" a 

where Jaa' — Ja'a initial exchange couplings in the interaction part Hex = X^qq' '^aa' ^aa' ^ Hamilto- 

nian, where s„„, = iE^j^^^,^^, a|;,^j^CT^^, OavA;'. 

The full RG equations at finite voltages and temperature are quite involved and will be discussed in detail in a 
forthcoming work^. Here we consider the equations either for /io, = or for T = 0, which simplifies them considerably. 
Furthermore, we have taken the special case of two reservoirs with /i^ = —fiR = and used the special initial 
condition Jaa' — '^\/XaXa' J with < < 1 and xl +Xfi = 1 (called "2-reservoir case" in the following). In this case 

the conductance in units of Go = ^ follows from G = G{E = 0) with G{E) = 7r^Fi(£'). In the general case, the 
variation of the current flowing in reservoir 7 is given in Laplace space by 5I^{E) = '^^tt5T^{E). The spin dynamics 
follows from {S){E) = -g^^(5)(t = 0). 

For zero temperature T = but arbitrary nonzero /ia a lengthy but straightforward calculation gives the following 
RG equations 

o2 c\ 

T{E) = zHi2Ji2(-B)J2i(^i2) , ^5T{E) - ^(5/il2Hl2 Ji2(-B) J21 (^12), (16) 



— {nSr^iE)) - -^Hi2/72(i^)i^2i(^)('5/ii2 + i^i2<5ri2), (17) 
d 1 1 

— Jl2(-E) = --H13 Ji3(i?) J32(£'i3) — -H32 Ji3(£'32) J32(£') 

— - Jl2(^^34)>^34(^')>^43(-£'1234)(-^34ni234 + ^12341134), (18) 
^^12(7;) = -H23Jl3(-B2i)if32(£^)-n3ii^i3(£;)J32(i;2l) 



dE 



+ ■^Jl2{E2134)J34:{E2l)K43{E){Z34ll2134 + -^21341134), (19) 



^ -Hi3/73(^^)^32(i?13)-n32Jl3(i?32)/3^2(^^) 



+ 2'^12(^34)^34(-E) J43(i?i234)(^34Hi234 + ^1234H34), (20) 
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where Z^i = Z[E^2). Z{E) 



-, STi2 = ST{Ei2), = n(£;i2) and ni234 = n(£;i234), with n(s) = e+^- 



These equations are even vahd for an arbitrary number of reservoirs and arbitrary initial matrix Jaa' ■ We emphasize 
again that in these equations the index 1 = ai contains only the reservoir label, and the new notations E12 — E + fj,i2, 
/^i2 = Ml — M2j and -E1234 — E + ^12 + ^34 have been introduced. 

For zero voltages = we consider the "2-reservoir case" and parameterize the renormalized vertices Ji2{E) 



2y/xT¥^J{E), KuiE) = 2y^^K{E) and I^2iE) = {6a^^ 
of all quantities, we obtain the following RG equations 



■'C(27 



)2y/xiX2Ji{E). Omitting the common £^-argument 



dK 
'dE 
dG 
'dE 



4z J2(^^(i) - 4JF(2)), 



d£ 
'dE 
dJi 
'dE 



_2J2(f(3)+ZJF(i)-2JF(4)), 



where we have defined the integrals 



F^^\E) 

f^^He) 
f(^\e) 



Z{E) du du 



Z{E) J duj I duj 
Z\E) 



2TTiT 

-Z{E) 



[1- 
duj 



uj + Lo' + x{E) 

,F {E,u:)nuj)nuj' 
uj + uj' + x{E) 
d^ 



Z{E) E 



2{2mTf 



{iuJn + iuJm + x{E)y^ 



Z{E) d^ 
2TTiT'dp' 



1 d^ 



Z{E) / duj 



Lo + x{E) 

F{E,Lo)r{Lo) 
uj + x{E) 



= z{E) y: 



2-KiT 



Z{E) d 
2TiiT dp 



ZHE) 
2mT 



iiln + 1 + p) 



and 



F{E,uj) = Z{E) J dw' riu') 

00 

= -Z{E) 2TriTY^ 



1 



uj + uj' + x{E) uj' - 
1 



X{E) 



with = fiuj) 



UJ + iuJn + x{E) iuJn + x{E) 
and the Matsubara frequencies w„ = 2'KT(n -\ 



(21) 
(22) 
(23) 



[P^{P)] , (24) 

(25) 
(26) 

(27) 



(28) 



As seen above, all these 



integrals can be expressed through the digamma function tp of the argument p — 



2iriT ■ 



In the following sections we will analyze these RG equations analytically in certain regimes. Furthermore, we will 
explain how to set up the initial conditions in order to enable an achievement of the scaling limit. 



A. T=V=0 



At T = y we use the RG equations (11ID-(I121) together with F^^) = i^^^) ^ and F^^) = F^'') = 0. With 

E ^ iA and A = A + r(«A) we obtain 



2r 2JH1 + ZJ) AKJ{1 ~ ZJ) 2JiJ Stt^ JiK 
Ar = — caJ = , dAK = , dAJi = —, dt^G = -A^XlXr—;^ 



(29) 



where Z — ^_^^^p . With J — ZJ, J/ — ZJi and I = In 4^ the equations can be written as 



'^ = -'^Zj\ ^ = 2P(1-J), ^^AKJ{l-Jl ^ = 2J,J(1-2J), '^^Axr^xn^-^JiK. (30) 
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These equations have the following invariants which are fixed as follows 

rr X 1 ^ o K Jj 

iK = ^e. 2j = — — 1 = = — 2 — —, 1 = ^= — . (31) 

(l-J)l/2 (l~J? J2 ' J{l~J) 

Hereby Tx sets the energy scale and defines the Kondo temperature. The other invariants are fixed by comparison 
with the initial conditions ([Tsl) in the scaling limit A — ?> oo and J — such that Tk is kept constant. The initial 
conditions (|15p are given hy Z ~ 1, K = 2J^, and J/ = J, which, in the scaling limit, lead to the invariants (|3ip . 
Using the results for Jj and K, the RG equation for the conductance can be integrated to 

G = AxLXB — J^ , (32) 

where the integration constant has been fixed again by comparison to the initial condition ([T5l) . Having fixed the 
invariants all information about the microscopic details of the model parameters is swept into the single energy scale 
Tk. The only remaining differential equation to be solved numerically \s d\\ = = ^j^^ , where J is a function of 

A via the expression for Tj^. To solve it we need the initial value A(0) = r(0) = F or, equivalently, the initial value 
J(0) = J. This value is determined from the condition that the normalized conductance is unitary G(0)/{4:XlXr) = 1 
a.t T — V — for the 1-channel Kondo model. This yields the values J = -^^^ 0.37 and T « 5.11 Tk- ^From the 
numerical solution we obtain all quantities at some high energy scale Aq ^ Tk- At finite T or /i^, we take these 
values as initial condition at Aq 3> T, \na\,TK and solve the RG equations starting from Aq. 

In second order truncation we neglect the terms ^ J^,KJ^ on the r.h.s. of and obtain Tk — XJe^^j^ 
Z = J ' ^ ~ 23"^ , Ji — J and G — Ax^xji^^j'^ . Fixing unitary conductance yields J = J(0) = w 0.37 and 
f = F(0) = w 10.57 Ti<- for the initial values for F and J at A = 0. 

Weak coupling- — To obtain the spin relaxation rate F(iA) in the weak coupling regime A ^ T^, we expand J 
and F around the 2-loop poor man scaling solution Jp(A), defined by In ^ = — h ^\tl^—^- Using F ~ AO{Jp) 

and J ^ Jp + 0{J^) we find A A(l + 0(Jp)) and i In ^ - i In = ^ + 0{Jp)- This can only be 
fulfilled if J = Jp + 0{Jp), from which we obtain d^X = ^ = = 1 + 2Jp + 3Jp + 0{Jp)- Using the ansatz A = 

A(l + aJp + 6j2 In Jp + cj2 + 0(J^)) and d^Jp = -2J2(1- Jp)/A, we find ^aA = l + aJp + 6j2 In Jp + (c-2a) J2 + 0(J3). 
Comparing the two expressions yields a = 2, 6 = and c = 2a + 3 = 7, i.e. with E = iA 



J{E) = Jp{-iE) + 0{Jl), T{E) = -2iEyp{-iE) + ^^JH-iE) + 0(j3)| . (33) 

Time evolution. — The spin dynamics is calculated from the inverse Laplace transform {S){t) = ^ / dES^{S){0), 
with X{E) = —iE + T{E), where the integration contour lies slightly above the real axis. The resolvent has a pole 
at A(£;) = 0, which corresponds to the fixed point J{E) = J* = I since A = Tft:e"37(l - J)i/2/ji/2_ ^he fixed 
point is reached for —iE — A* — — F* < since at A = we have chosen J(0) — — ^ < 1. Close to the fixed point 
A - A* = -F* we get ^ = AZp/X « aZ/X with a = 4J*^ This leads to Z - A". Using ^ = |- A"" we get 
A - (A- A*) — . This means that the resolvent has a branching pole at i? = —iV* followed by a branch cut on 
the line E = — z(F* + x) with < a; < oo. The contour integral can be closed in the lower half around the branch cut 
and, in the long-time limit t ^ the spin dynamics follows from 

^ -r** / j„„-2:tT /c\/n^ _ „-r*t / „-xt ^ _ V „-r*t 



with the exponent 



4j*2 



1 + a l + U 



.2 • 



(35) 



As a result the exponent depends on the fixed point value J*, which in turn depends on the truncation order. In 
third order truncation we get J* = 1 leading to g = |. In second order truncation we neglect the term of 0{J^) 

in the RG equation ([29|) for J. This yields ^ = 2J^(1 — 2J) and Z = 1 — 2J. The fixed point in second order 
truncation is at J = J* = ^ which yields g = i. A fixed point value J* — oo would lead to g = 1. This can be 
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achieved by changing the /3-function by hand to ^ = j— , which is still consistent with the /?- function in third order 

truncation. This leads to the invariant Tr- = AJ^/^e~2j and to Z = e^'^'^^'^ . Close to the fixed point A = 0, the 
coupling J ^ jj- tends to infinity and fulfills the differential equation d\J ~ ^'k\'\^ ~ _j3'/2^2j+j _ jj-^ ^^^^ leading 
order we get J '-^ (— In Y^"^ and A ~ J^^/^ ^ {— In )^^^^, which results in a pre-exponential function 
~ 7(h7T^l)374- Thus, besides the power- law j the pre-exponential function contains an additional logarithmic term 
in this case. In conclusion, the pre-exponential modulation depends crucially on the form of the /3-function and can 
not be determined within a perturbative truncation scheme. Instead, one needs a non-perturbative analysis of the 
/3-function close to the fixed point which goes beyond the scope of this work. 

The problem of obtaining the correct pre-exponential modulation for the 1-channel Kondo model is analogous to 
the problem of obtaining the correct scaling dimension close to the fixed point for multi-channel Kondo models. For 
channel number > 1 the RG equations for Z and J change to ^ — —4NZJ'^ and ^ — 2J^(1 — NJ), giving 
Tk = Ae~2j( ^ Afj )^^^ ^^"^ ^ = (1 ^ NJY. The fixed point A = corresponds to J = -i- and the system should 

reveal non-Fermi liquid behavior at low energies^. This means that we have to choose J{E = 0) = identical to the 
fixed point value such that the pole A = lies at the origin £^ = A* = of the complex plane with r(0) = F* = 0. In 
contrast to the exponential decay for the 1-channel case, this leads to a pure power-law decay in the long-time limit, 
which is another signature of non- Fermi liquid behavior. Expanding around the fixed point we get ^ ~ aZ / A with 

a = ANJ*'^. According to (l35t this gives the power-law exponent g = = j^^^jn- Inserting J* = we obtain 
g — However, this result is only consistent for ^ 1, where the /3-function can be systematically truncated 

and the system stays in the weak-coupling limit^. For small N the result is not reliable like in the 1-channel case 
discussed above. 



B. T=0 and finite voltage 



At zero temperature we start from the RG equations (|16II20[) and discuss analytically the "2-reservoir case" in the 
weak coupling regime of large voltage V » Tk and the strong coupling regime of small voltage V <^ Tk- 

Weak coupling. — To obtain the weak coupling expansion for V 3> Tfc we follow Ref. |3 and discuss the solution of 
the RG equations separately for A > V and A < V, where A = —iE. For A > T^, we expand the solution of the RG 
equations around a reference solution, where all propagators are replaced by ^. No voltage dependence is generated 
in this case and the RG equations obtain the same form as in ([50)1 but with / — 7> In As a consequence, the reference 
solution is given by pi|) but with J — > Jp, where Jp{A) is the 2-loop poor man scaling solution, already introduced 
in section Hi A1 To calculate the first correction to the reference solution for A > 1^, we insert the exact identity 

-^^A = -itHa + — < Zi2(ln )iJA 



Ai2 + ri2 A dAi A 

- (In ^^^^)-^ {ZuH,} + ^^H, (36) 

for all propagators occurring in the RG equations (jl6ll20[) . where A12 = A — 1^12 and 77a stands symbolically for 
the rest of some term on the r.h.s. of the RG equation involving products of various couplings. We get Hf^ ^ 
{H\ ^ Jp) for d{^J,dp^n and dj^Z {dj^K and dp^G, with G = -k-^Tl). The first term of (p6|) generates the reference 
solution, which, according to (I5T]) and ([5^ . gives 0{Jp) {0{Jp)) for J, H and Z — 1 {K and G). We are aiming 
at calculating all quantities one order beyond the reference solution. To estimate the various terms of (j36p we use 
Ai2+ri2 ^ Q(V_^ r ) _ Q^v_^ fQj. A > F, F, where we used F - 0{AJp) according to As a consequence, the 

third term of ([211) is - Oi^J^, \J^) 0{^J^, i J^)) for d^J^dAl^ and d^Z {d^K and d^G). Using ff ^ j^Oi^), 
the fourth term of (|36| is ~ 0{^Jp) J^)). Thus, after integrating over A, the third and fourth terms of (|36)) 

produce corrections ~ Jp Jp) for J, P and Z {K and G), i.e. two orders higher than the reference solution. In 
contrast, the second term of ([55]) produces the first correction ^ ^O(Jp) (~ \0{Jp)) to the reference solution. We 
can neglect the influence of this correction on the first term of (|36p . since this gives contributions of the same order 
as the third term of Thus for our purpose it is sufficient to take into account the first two terms of and 

calculate Hp^ by using the reference solution. For the RG equation ((T7)) for the conductance this means that the term 
involving (5Fi2 ~ SVJp can be neglected on the r.h.s. and the first term on the r.h.s. leads to the first correction to 
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the reference solution, i.e. the weak coupling expansion in the regime A> V reads 

GiE = V) = AxLXn^J^iA) { 1 - 2Jp(A) > ; In lli^-Ilif + 0{J') } . (37) 



Following Ref. we use this result at A = as initial condition to solve the RG equation for the conductance in the 
following regime A < y by expanding in the reference solution Jp{V) evaluated a.t A = V. Replacing all vertices on 
the r.h.s. of the RG equation for the conductance by the reference solution a.t A = V, we get with Hy ^ JiK ^ 0{Jp) 
and if ~ A^O{J^) instead of m the expression j^^^Hy = ^ {Zi2(ln Aia+Tia + ^0(j5). Applying 

Jy dA to this expression and adding the initial condition ([57)1 evaluated at A — V, the second term on the r.h.s. of 
(137)) obviously cancels against the contribution from the lower bound A = y of the integration. The remaining first 
term of ([37l) evaluated a.t A = V together with the contribution from the upper bound A = of the integration gives 

for the conductance at E = the result G = Axlxr^J^{V) {l - 2Jp{V) Eai^a^ In '"^"^'^^^^'^ + where 

we used ^i2|a=o = 1 + 0{Jp{V)). With /ii2 ~ ai2-j and ai2 — ai — a2, we find 

E 1- V"^^^"^ - E ln(-4"-) + O(^) = In(-z) + ln(^) + 0(1) = 0(I) ^ (38) 

where we have used T(E = 0) ^ 0{VJp) which can be derived from an analogous weak coupling expansion. Thus, 
we can neglect this term and finally get the weak coupling result 

G{V) = AxLXR — J^iV) + 0{4). (39) 

Strong coupling. — In the strong coupling regime V <^ Tk the 1-channel Kondo model behaves like a local Fermi 
liquid^, which, in particular, implies a quadratic dependence of the differential conductance on the applied bias 
voltageSiia G = 4:Xlxr{1 - cviV/T)^ + 0{{V/r f)}. As a reference scale we use here f r{E = 0). In the following 
we recover it by means of the RTRG as well as provide a quantitative value for the Fermi-liquid coefficient cy up 
to the next-to- leading order in the renormalized coupling constant J = J{E — 0)\t=v=o- This expansion can be 
performed by differentiating the RG equation (fTT]) for the normalized conductance twice by the voltage and integrate 
it over E. A more direct path is to consider the original diagrammatic series and calculate directly the third order 
variation (5^E-^(_E) = ■^-^yj'Ej{E)dV^ of the current kernel by using the diagrammatic method explained in SectionH] 
Of particular advantage is the fact that the variation ^r(£') of the Liouvillian vanishes at y = T = for any E. 
This follows directly from the RG equation (IT6l) after inserting the algebra of the vertices at V = T = derived in 
Section Hi Al given by 

Jaia^ 2^/x^aiXa2 J ^ ^aia2 2 , ^aia2 — {^aij ^ck27 )^'\/*^cki *^ck2 

Z={l-Jf, K = 2P, Ji = J{l-J), (40) 

where J = ZJ and J/ = ZJj. As a consequence we can express the variation of the propagator at y = by derivatives 
with respect to the Laplace variable S"Ili2 — -^{SfLi2)"d^Ili2- With this property we get for the third order variation 
of the current kernel at zero voltage 



2 3! 

1 2 

+ hsp,2rsfi^, i^-z^rw^r- H + bfi^2{Sfii3)' i^^^^l + o{j% m 

^ 123 ^ 123 

where all diagrams have to be evaluated at y = on the r.h.s. The corrections are of 0{J^) since for the current 
kernel one of the vertices has to be ii' ^ J^. At y = the last two diagrams are identical and, by using partial 
integration, can be written as 

12 3 12 3 ^123 
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Analogously to ([5]) we can perform two partial integrations with the first term on the r.h.s. of (|4T]) and, by again 
using P^ . we arrive at 

a3 1 I ^ 1 



1 2 



1 



12 3 

Due to the derivative of the Fermi functions in all contractions, all frequencies are set to zero and the diagrams can 
be straightforwardly evaluated by performing the sums over the indices rj and a by using the algebra of the vertices 
together with their value (HH)) at y = 0. With fia — \oiV , ai2 — ax — cl2 and 8ax-^ ~ ^a^-^ = 570^12, we obtain after a 
straightforward calculation 

xJ^(ai2)^+ ai2((ai2)' + (ai3)'-3(ai2fai3-3ai2(ai3)') i'J'^' + O(J'). (44) 

Performing the sums over and using 11^{E) = iT^{E)L'' along with G = tt^F^ and (|40l) . we obtain at = the 
final result 

Finally, we note that the second variation of the current kernel at = yields zero since an analogous derivation 
as above gives 

12 12 3 

= i(5/li2)' - ^(5u„-Su,,)' i^^j^i + (46) 

#9 •• 1 2 3 



Inserting the algebra of the vertices along with their values (|40l) at = 0, and using X]aiQ2 ('^i^)^ — and 
X^aiaac, ("-"^is ~ cn3)^cn2 = 0, wc find & vanishing value. This shows that the corrections of the conductance to 
the unitary value start with a term ^ V'^ as required by the Fermi liquid theory. 



C. V=0 and finite temperature 

At zero voltage we consider the " 2- reservoir case" and start from the RG equations (|21j|23p to discuss analytically 
the weak and strong coupling regimes T » Tk and T ^ Tk- 

Weak coupling. — To obtain the weak coupling expansion for T ^ Tk we proceed analogously to the derivation 
presented in Section Hi Bl We set E — iA and start with the regime A > T, where we expand in the reference solution 
Jp(A). As already discussed in Scctions llI Bl and lll Al the reference solution pip with J Jp is obtained if one replaces 
the integrals F(')(iA) defined in (|24j|27p by iF^^'^) | and F^^"* ^ q in the RG equations (|21j|23l) . Therefore, we 
construct an analog of §^ by writing iF'-^'^^ = ^ + ZdAp'-^^^^ +0{J^/A), such that F^^'^) _ 0{{^ f, Jp) for A > T. 
Using p — ^2^y^ and |^ — 2^i^ + 0( J^)) together with and we find that this can be achieved by defining 
F(i) = ^{p5f(p)} _ 1 _ In ^ and F^^) = ^-(p + i) _ In Furthermore, we write iF^^.^) = ZcLaF^'^-''^^ + 0{J^/A) 
such that F(^^^) ^ 0(^)^, which follows directly from the asymptotic expansion 

F(2.4) ~ ^ for a > T. The analog 

of ([M]) then reads 

zF«i/A = c^i/A + ^ {zAFf I/a} - Ff ^ {ZaI/a} + 0(1 J^i?), (47) 

with ci = C3 = 1 and C2 = C4 = 0. All arguments used in Section [II BI can then be overtaken in the same way, one just 
has to replace the logarithmic function In ^i^+^i^ by the functions F^*' in the corresponding terms. In particular this 
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means that the terms with the integrals F^^'"^^ do not contribute to the first correction beyond the reference solution 
since they appear in third order 
by analogy with ^ for A > T 

^(A){l-4J,(A)F«(r) + 0(j2)}. (48) 



since they appear in third order and fall off sufficiently fast ^ for A ^ T. Thus, for the conductance, we obtain 

37^2 

G{E = tA,T) = 4xLXR — J^, 

Following Section lll Bl in the next regime A < T, we expand in the reference solution Jp{T) evaluated at T. Replacing 
all vertices on the r.h.s. of ([25]) by the reference solution at A = T, we see that only the first term involving F'-^'> can 
correct the result in 0{Jp). Furthermore, for the integration Jj, dA, we use the identity iF|^^ = ZA^rdApi^^ +0{^J^), 



with fI ' = ^{/9*(/5)} - 1 - In ^21^ which can be proven by using \Za - Za=t\ 0{J^), 



dA 



^(l + 0(j2))and 



F^^'' ^ 0{^) for A < T. This form is constructed in such a way that pj^^ is well-behaved for A ^ and, at the lower 

boundary A = T of the integration, is exactly identical to -^^=7-' such that the third order term of (|48l) is canceled. 
The remaining first term of (^5)) evaluated at A = T together with the contribution from the upper bound A = of 



the integration gives for the conductance at E' = the result G = AxLXR^^Jp{T) |l — 4Jp(r)/'|^p + 0( Jp)|. Using 



rA=o ^ 0{TJp) (which can be derived from an analogous weak coupling expansion) we get pa=o — ^'^^"Jt^" 0{Jp), 
i.e. A = is equivalent to p = up to higher-order terms. Using j^{p'^{p)}p=o — —J, where 7 = 0.577 ... is Euler's 
constant, and Za=t = 1 + 0{Jp), we get -Pa=o ^ ^7 ^ 1 + In (27r) + 0{Jp), leading to the final weak coupling result 

G(r) = 4xLXR^J^iT){l - Jj,(r)4(ln(27r) 



1 



7) + 0(4)}. 



(49) 



Strong coupling. — In the strong coupling regime T <^ Tk we aim at deriving the Fermi liquid coefficient ct in 
the low temperature expansion of the conductance G — AxlXr{1 — ct{T/T)'^ + 0{{T/T)^)} or of its first derivative 
= —Sxi^xryjCtT + 0{T^ /T^). To determine the latter quantity we consider the T-derivative of the RG equations 



2TII23p and integrate them formally over E (alternatively, one can also derive a diagrammatic technique for T- 
derivatives similar to the one described in Section U for E- and voltage derivatives^) . Employing the representation 
F^^){E) = ^^ij(*)(2^) we find the useful relations 



ai^(i/3) 
dT 



Z_d_ 

TdE 



^i^(l/3) 



z 



0(J2 



dT 



Z_d_ 

TdE 



^i^(2/4) 

~~z 



0(J2 



(50) 



along with the asymptotic behavior xF^^'^^ — Z , ;^E'(2,4) ^ O(^) for E :$> T. More precisely, for E ^ T, the 
corrections to (|50|) are of 0{J'^ / E"^), such that they can not lead to a reduction of the order in J when integrating 
the RG equations in the regime from E = 00 to E = T. This follows from analyzing the result of the integration 



procedure, which shows that, for i? >• T, we get "^'^q^'^' ^ -^J^ ■, ^ W"^^^ ^nd |y 



d(J.,J,.Z) T t2 dK 
8T ^ ' dT 

made use of ^ = 1 + 0{J^). As a consequence, one finds for the T-dcrivative of ([22] 

d d{JiK) 6J1KJZ d fxF^^^ 



1 



^J^. Furthermore, we 



(51) 



dE dT T dE \ Z 

where the corrections are again of 0{J^ / E"^) for E ^ T. Integrating ([51]) over E we obtain to the leading order in J 

d{JiK) _ 6J1KJ 



- Z 



o{.r 



dT T 

Using this result along with ([50[) to calculate the T-derivative of (l23l) we obtain 



d dG 
dEdf 



Gtt^xlXr, 
T 

en'^XLXR 



d_ 

'dE 



T 



d_ 

'dE 



JiK{xF^^^ - Z)- GJjKJxF'-^'^ 



z 

6JIKJZ 



1 - 6 J 



Z 



where we used ^^'gj^^ 



dE 



QJiKJF^^^ -f 0(J5) and |f - 0{J^) in the last step. Using 
;e this result over _ 

dG Gtt^xlXr 



+ 0{J^), we can integrate this result over E and obtain the equation 



dT 



T 



-JiK 



{xF^'' 



Z] -6J 



{xF 



(2) _ ^(2) 



OiJ% 



(52) 
+ 0(J') 

+ 0(J5) , (53) 
and ([^B| as well as 

(54) 
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where ^ ^2 



1 



with the asymptotic behavior F^'^'^ '-^ 0{^) for E ^T. 



Setting £' = in ([M)) and using the low temperature expansions x^'^^^ ~^ ^(1 ~ 6p^)' X-?"'^'' 
24p2 p - 2^ > 1, we get 

dG Att'^xlxrJiK 



72 1 



(Zf)2 



[1 + 3J) T + 0{T^,.r 



and 



(55) 



where Z ^ Z(0), T = r(0), J = Z(0) J(0), J/ = Z(0) J/(0) and X = K{Q). This expression leads us to the final resuh 
for the Fermi liquid coefficient ct 



ct 



2 Z2 



(1- J)3 



(56) 



where we used J/ = J(l — J), i^T = 2J^ and Z = (1 — J)^ [cf. (|3ip ] in the last step. Comparing ct with the Fermi 



liquid coefficient cy [cf. pS)) ] we get the familiar universal ratio 
Kendo model. 
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